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Hard diffraction in events witli dijets and rapidity gaps has been studied by D0 and CDF 
for three processes: hard color singlet exchange, hard single diffraction, and hard double pomeron 
exchange, using Tevatron pp data at — 630 GeV and 1.8 TeV. Measurements of rates, rj, Et and 
y/s dependencies are presented and comparisons made with predictions of several models. 
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I. INTRODUCTION 



Events with a region of rapidity space devoid of particles (rapidity gaps) were first observed in cosmic ray data 

0. The idea of diffractive dissociation of projectile and target to produce such events soon followed Later the 

, interpretation of total, elastic and diffractive cross sections in terms of the exchange of an object with the quantum 

' numbers of the vacuum, called the pomeron, proved very useful ^ . F. Low and S. Nussinov suggested that the pomeron 

corresponds to the interchange of two gluons M. In a latter development Ingclman and Schlein (IS) proposed that 
T~H ' . . . 

^ ■ high pt jets could be diffractively produced via pomeron exchange and that this might probe the partonic structure 
^\ . of the pomeron . Events containing rapidity gaps and jets were first observed by UA8 |^ , giving rise to the field 



of hard diffraction. This area of interest has expanded considerably in the last decade, with the availability of high 
energy colliding beams. Jet production with rapidity gaps have been observed at the Tevatron and at HERA 
[H . CDF and D0 have studied dijet + rapidity gap events using data from the 1992-1996 Tevatron pp collider run 
(Run I) at center-of-mass energy (-^/s ) = 1.8 TeV, with a short run at 630 GeV during that period. Because of the 
5 ' large center of mass energy and large integrated luminosity, the new CDF and D0 results can give further insight 
into diffractive processes. 



The D0 detector is described elsewhere |l(l. Jets are found in the uranium- liquid argon calorimeters using a cone 



O . algorithm with radius R — 0.7 in the rj ~ (f) plane ]11[ |. Particle multiplicity is determined in the central region 
' {\rj\ < 1.0) using the number of towers (0.1 x 0.1 in A77 x A0) with transverse energy {Et) above 200 MeV in the 
. ^ ' central electromagnetic calorimeter and the number of tracks in the central drift chamber. In the forward region this 
multiplicity is measured by the number of towers with (Et) above 125 MeV in the electromagnetic end cap calorimeter 
(2.1 < I77I < 4.1) and 500 MeV in the hadronic end cap calorimeter (3.2 < \t]\ < 5.2). Because the last layer of the 
hadronic calorimeter (at the limit of the forward acceptance) is composed of stainless steel and produces less noise 
than the uranium sections, the threshold for particle detection was reduced to 50 MeV in this layer. In addition, in 
the forward region we also use an array of scintillator hodoscopes, called L0 detector, to tag the presence of charged 
particles in the region 2.3 < < 4.3. 

The CDF detector, described in ref. |lj], consists of a large central detector with tracking in a solenoidal field and 
calorimetry over jryj < 4.2. To measure particle multiplicities, CDF uses the central tracker {\ri\ < 1.1, pj'^'^"-'^'^ > 300 
MeV), the central calorimeter {\ri\ < 1.1, Et*'°''"'^^ > 300 MeV corrected), and forward calorimeters (2.2 < |?7| < 4.2). 
For the last two months of the collider run, CDF installed three Roman Pot detectors to trigger on quasi-elastically 
scatted antiprotons. Nearly all the pot triggers have 0.05 < ^ < 0.1, where ^ = 1 — ccf is the fraction of momentum 
lost by the antiproton and carried by the pomeron. 
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II. HARD COLOR SINGLET EXCHANGE 



Two jets separated by a rapidity gap has been proposed as the signature of color singlet exchange (CSE) carrying a 
high [ p^JTil . Rapidity gaps between jets have been observed both at the Tevatron Q and at the DESY ep Collider 
(HERA) Q. The measured rates of w 1% at the Tevatron and w 10% at HERA are too large to be accounted for by 
electroweak boson exchange and indicate a strong interaction process. 

D0 and CDF have made recent studies of dijet data with central rapidity gaps. Both experiments measure the color 
singlet fraction (/g) at of 630 GeV and 1.8 TeV. The observed color singlet fraction includes the probability that 
the rapidity gap is not contaminated by particles from spectator interactions. This survival probability (S ^ 10% at 
1.8 TeV) is assumed to be independent of Bjorken x and the flavor of the initial partons in the hard scattering 18 
but depends on 4~s (S'eso/S'isoo = 2.2 ± 0.2) ||l|l. 

CDF measures the fraction of colorless exchange to all opposite side dijets from the tracking distribution. Results 
are listed in Table ||. The ratio of the CDF fractions from the measurements at the two center of mass energies is 
_R(r^^) = 2.4 ± 0.9. No Et dependence of the signal is observed. 




FIG. 1. The multiplicity between the dijets for the D0 high-i^T 1800 GeV sample: (a) two-dimensional multiplicity, n^a.i 
vs. ntrk\ (b) ricai only with NBD fit, plotted on a log-log scale to emphasize low multiplicity bins. 



experiment -y/s (GeV) jet Et t?"''^' nb. of triggered events 

CDF 630 > 8 GeV 1.8< |r?| <3.5 Ik Ik(same-side) 

D0 630 >12 GeV I77I >1.9 7k 

CDF 1800 >20 GeV 1.8< |j7| <3.5 10k +30k(same-side) 

D0 1800 >12 GeV |7?| >1.9 48k 

D0 1800 >25 GeV |r;| >1.9 21k 

D0 1800 >30 GeV H >1.9 72k 



TABLE I. Kinematic cuts for color singlet exchange. 



./s(% ) 


Et (GeV) 




experiment 


2.7 ±0.7 (stat.) ±0.6 (syst.) 


> 8 


630 


CDF 


1.85 ± 0.09(stat.) ± 0.37(syst.) 


> 12 


630 


D0 


0.54 ± 0.06(stat.) ± 0.16(syst.) 


> 12 


1800 


D0 


0.94 ± 0.04(stat.) ± 0.12(syst.) 


> 30 


1800 


D0 


1.13 ± 0.12(stat.) ± 0.11(syst.) 


>20 


1800 


CDF 



TABLE II. Color singlet Fractions . 
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D0 has recently published results for dijet + central gap events (See Table | for kinematic cuts.) Single 

interaction events are required. The particle multiplicity in the central rapidity region is approximated by the 
multiplicity, ricai, of transverse energy above 200 MeV in the electromagnetic calorimeter, and by the track multiplicity 
in the central tracking chamber, ntrk- Figure 1 (a) shows the D0 multiplicity distribution for ricai versus ritrk- 

To calculate the fraction due to color singlet exchange, the leading edge of each ricai distribution is fitted using a 
single negative binomial distribution (NBD). The fraction of rapidity gap events (/g) is calculated from the excess of 
events over the fit in the first two bins {ricai = or 1) divided by the total number of entries. Figure 1 (b) shows the 
ricai distribution and the NBD fit for the high Et sample. See Table |l| for values of fs- The D0 value of the ratio 
of the rapidity gap fractions at 630 and 1800 GeV is R{j^) = 3.4 ± 1.2. 

Measuring the color-singlet fraction as a function of Et, f] and y/s probes the nature of the color-singlet exchange 
and its coupling to quarks and gluons. If the color-singlet dynamics are similar to single gluon exchange except for 
different coupling factors to quarks and gluons, the color-singlet fraction would depend only on parton distribution 
functions via xp- Thus for a color-singlet that couples more strongly to gluons than quarks, the color-singlet fraction 
would fall as a function of increasing x, since the gluon distribution becomes suppressed relative to the quark distri- 
bution as X increases. This implies a decreasing color singlet fraction with increasing jet Et and A77 or decreasing 

To measure the color-singlet fraction as a function of Et and 77, D0 uses the two-dimensional multiplicity, 
(jicai vs. ritrk) which gives improved signal-to-background ratios compared to the NBD method. This is useful 
for smaller statistics samples and avoids large uncertainties in the color-exchange background subtraction. The "2D" 
color-singlet fraction /2D is defined as the fraction of events with ricai + ntrk < 2. The results are shown in Figure ||. 
The systematic errors include effects from background estimation. The measured color-singlet fraction shows a slight 
rise as a function of dijet Et and Arj. 

To compare the experimental color singlet fractions to models, D0 uses HERWIG 5.9 which includes a two- 
gluon exchange with BFKL dynamics jl5|, and uses CTEQ2M parton distribution functions. In addition, D0 uses 
the t-channcl photon exchange process in HERWIG to investigate models in which the color singlet couples only to 
quarks with a massless photon-like singlet. 

In the soft-color rearrangement model p^], initial state quarks have fewer color combinations and thus, a higher 
probability of being rearranged into a colorless state, than initial state gluons, i.e. Cgg < Cqg < Cqq, where the "Cq;,s" 
are the effective color factors representing the couplings to different initial state partons. A reasonable choice of color 
factors is Cqq = i, Cqg — ^ and Cgg = (^). Predictions of these models are simultaneously fit to the experimental 
Et and A77 dependence of fs at ^/s— 1.8 TeV, letting the normalization float. The results are shown in Figure ||. 
The data favour color-singlet models that couple more strongly to quarks than gluons, but a single- gluon model (no 
dependence) can not be excluded. 
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FIG. 2. Monte Carlo fits to the measured color-singlet fraction /2D. The normalization is allowed to float. 
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III. HARD SINGLE DIFFRACTION 



In the IS picture of hard smgle diffraction, a pomeron (color singlet object) is emitted from the incident p (p) and 
undergoes a hard scattering with the p (p), leaving a rapidity gap in the direction of the parent particle. The signature 
is two jets produced on the same side and a forward rapidity gap along the direction of one of initial beam particles. 
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CDF PRELIMINARY 
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FIG. 3. Ratio of Data to Monte Carlo simulation as a function of /3 using a flat gluon distribution and standard pomeron 
flux 



Data was taken by the CDF detector at the end of Run I using a trigger which requires tagging the recoil p with 
"Roman Pot" detectors. The typical acceptance for these detectors is 0.05 < ^ < 0.1 and < |t| < 2 GeV^, where 
^ is the fractional momentum lost by the antiproton and t its four momentum squared. After applying several cuts 
to select events with a good reconstructed track in the Roman Pots, CDF extracted the momentum fraction of the 
interacting parton in the pomeron, /3, for dijet events with Et > 7GeV, using the following expression: 



= 



be am 



(1) 



The (3 distribution for the pomeron was obtained by subtracting several background contributions in the data, of 
which the most importants are 1) non-diffractive dijet events accidentally overlapped with a Roman Pot hit, 2) meson 
exchange background and 3) double diffraction background. After subtracting these contributions from the data, then 
unfolding the detector acceptance by using simulations with a flat gluon distribution, the data was divided by Monte 
Carlo simulations based on POMPYT |22j using a flat gluon distribution and the standard Donnachie and LandshofF 
flux parametrization [2^ . The comparison, shown in Figure ^, shows agreement in shape for [3 > 0.2, but there is a 
discrepancy in the normalization by about a factor of 6, as well as an enhancement for the low /3 region. 
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1 800 Forward Jets 1 800 Forward Jets 




630 Forward Jets 630 Forward Jets 



FIG. 4. Number of hits level (iilo) and calorimeter towers (hcal) above threshold opposite to jets in the forward trigger 
sample at center of mass energies 1800 and 630 GeV. 

The D0 data were obtained using an inclusive jet trigger or a forward two jet trigger. Events are selected with two 
jets with Et > 12 GeV and \r]\ >1.6. The number of end cap calorimeter towers (ricai) above threshold is measured 
opposite to the leading two jets. The ricai distribution for the forward trigger sample at center of mass energies of 
1800 and 630 GeV is shown in Figure ^. A clear peak is seen in the ulq = ricai = (zero multiplicity) bin as 
expected for a diffractive signal. A two-dimensional fit on the tilo vs. ricai distribution, where data and background 
are fit simultaneously, allows the direct extraction of the fraction of events containing a rapidity gap at both energies. 
This gap fraction, including statistical and systematic uncertainties, is determined to be 0.64 ± Q.Qb{stat. + syst.)% 
for the 1800 GeV data and 1.23 + 0.10 - 0.09{stat. + syst.)% for the 630 GeV data. Work is in progress for the 
extraction of the pomcron parton distributions. 

IV. HARD DOUBLE POMERON EXCHANGE 

Central dijet events containing two rapidity gaps or a rapidity gap on the opposite side of a quasi-elastically scattered 
anti-proton have been studied by D0 and CDF respectively. This event topology is consistent with double pomeron 
exchange (DPE). The data can be used to give more information about the hypothesized pomeron. 

CDF took data at 1800 GeV center of mass energy using the Roman Pot trigger to tag antiprotons. A sample 
of 27,000 events with a tagged p and at least two jets with Et > 7 GeV is obtained. Low multiplicity events are 
selected by requiring NBBc{west) < 6 in the 16 element Beam-beam counters (BBC) on the same side as the pots. 
This gives 22,304 PJJ (pot-jet-jet) events. A sample of minimum bias events, with the same dijet selection, is used 
for comparison. Figure ^ (a) shows the calorimeter tower multiplicity in the east side (opposite to the pot track), 
NpcAL, versus the number of hits in the BBC. 
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CDF Preliminary 



CDF Preliminary 





^aj ''^^ °° (b) 

FIG. 5. (a) East side BBC vs Cal. tower multiplicity distribution for PJJ events, and 7 GeV dijets. (b) East side BBC vs 
Cal. tower multiplicity distribution for Monte Carlo Double Pomeron (with a flat /3 distribution.) 



In order to understand the shape of this distribution, Monte Carlo events were generated using a version of POMPYT 
modified to include double pomeron exchange, where the incoming proton and antiproton emit pomerons with a 
standard flux (Donnachie and Landshoff form [E3[ with parameters as measured by CDF 1 24 ) . The pomeron-pomeron 



interaction is treated like a hadron hadron collision which produces jets. A flat /3 distribution of partons inside the 
pomeron was assumed. Diffractive deep inelastic scattering data from HERA suggest such a hard structure with a 
rather flat /3-distribution |2^. The simulated events are shown in Figure || (b). The strong signal in the (0,0) bin 
only contains 24 % of the DPE events with £_ < 0.1 for 7 GeV dijets. 

By extrapolating linearly into the Nbbc — NpcAL = bin along the diagonal axis, CDF obtained the ratio of 
dijet gap to dijet no-gap events to be: 



i?( 



PJJG 
PJJ 



) = [0.36 ± Omistat.) ± {Om{syst. 



where G means gap (no detected particles for 2.4 < rj < 5.9), J J means two jets with Et > 7 GeV and P means a 
pot track with 0.05 < C < 0.1. When the discrepancy factor, D=18%, found in previous analyses of diffraction in pp 
and ep collisions 0, is applied (squared), the ratio from simulation is in good agreement with the data. 




FIG. 6. (a) Multiplicity distribution of calorimeter towers (ricai) opposite a tagged rapidity gap for the 630 GeV D0 data. 
The bottom plot shows an expanded view of the low multiplicity region, (b) Et distributions of the leading two jets for three 
data samples at 630 GeV. An inclusive sample requiring \rijet\ < 1.0 is shown in the solid histogram. The distribution with the 
added requirement of one forward rapidity gap is shown with dotted lines and the distribution for double gap events is shown 
in circles. 

D0 has taken inclusive jet data with a special trigger and searched for dijet events with two forward rapidity gaps 
along the direction of the proton and antiproton. Events were selected having two jets with Et > 12 GeV, \rijet | < 1.0 
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and a rapidity gap in the region 2.5 < I77I < 5.2. The multiphcity distribution of calorimeter towers and Level0 hits 
(71^0), on the opposite side to the rapidity gap, for data taken at 630 GeV center of mass energy is shown in Figure ^ 
(a). A clear peak at low multiplicity is observed above a fairly flat background in qualitative agreement with that 
expected for double pomcron exchange. Figure |^ (b) shows the Et spectra for the two leading jets in an inclusive 
sample with two central E-r > 12 GeV jets , a single forward gap sample and the double gap sample at 630 GeV. All 
three spectra are in good agreement where data are available, implying that the dynamics of leading jets produced 
in rapidity gap events appear similar to those of inclusive QCD production. Similar results are seen by D0 in data 
taken at 1800 GeV. 



V. CONCLUSIONS 

Recent studies of hard diffraction at the Tevatron have given new information about rates of diffraction and 
dependencies on Et, rj and y/s . 

The fraction of dijet events produced via hard color singlet exchange is about 1% at 1.8 TeV and is larger by a factor 
of 2 to 3 at y/s= 630 GeV. D0 has compared the Et and A77 dependence of the fraction of hard color singlet events 
to several models. The data favor a soft-color rearrangement model preferring initial quark states over two-gluon 
color-singlet models. 

CDF has preliminary results on the momentum distribution of partons in the pomcron using "Roman Pot" detectors 
to measure quasi-elastic scattered p in hard single diffractive events. D0 has studied hard single diffraction in forward 
dijet events and new results on pomcron parton distributions will soon be available. 

Both CDF and D0 have preliminary evidence for events with a hard double pomeron exchange topology. CDF 
has measured the fraction of pot dijet gap events to be 0.36% of the pot dijet events at ^/s — 1.8 TeV. D0 has 
studied gap-dijet-gap events at y/s= 630 GeV and 1.8 TeV. The Et distribution of the leading jets in double pomeron 
exchange type events is similar to other processes producing jets. 

More detailed studies of the Tevatron data and further comparisons with models should give more insight into the 
nature of hard diffraction. 



[1] P. Ciok et al, Nuovo Cimento 8, 166 (1958); K. Niu, Nuovo Cimento 10, 994 (1958); G. Cocconi, Phys. Rev. Ill, 1699 
(1958). 

[2] M. Good and R. G. Walker Phys. Rev. 120, 1857 (1960). 

[3] For a review see K. Goulianos, Phys. Reports 101, 169 (1983). 

[4] F. Low, Phys. Rev. D 12, 163 (1975); S. Nussinov, Phys. Rev. D 14, 246 (1976). 

[5] G. Ingelman and R Schlein Phys. Lett. B 152, 256 (1985). 

[6] A. Brandt, et al, Phys. Lett. B 297, 417 (1992). 

[7] D0 Collaboration, S. Abachi et al, Phys. Rev. Lett. 72, 2332 (1994). CDF Collaboration, Phys. Rev. Lett. 74, 855 (1995); 

D0 Collaboration, Phys. Rev. Lett. 76, 734 (1996); CDF Collaboration, Phys. Rev. Lett. 80, 1156 (1998). 
[8] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 79, 2636 (1997). 

[9] ZEUS Collaboration, M. Derrick et al., Phys. Lett. B 369, 55 (1996); ZEUS Collaboration, M. Derrick et al, accepted for 
publication in European Physics Journal (1998). 
[10] D0 Collaboration, S. Abachi et al, Nucl. lustrum. Methods A 338, 185 (1994). 

[11] The pseudorapidity, r] is defined as — lntan(S/2). The polar angle, 6 , is defined relative the the proton beam direction, (j> 

is the azimuthal angle, and Et is the energy measured transverse to the beam. 
[12] CDF Collaboration, F. Abe et al, Nucl lustrum. Methods A 271, 387 (1988). 

[13] Yu. L. Dokshitzer, V.A. Khoze and S.L Troian, Proceedings of the 6th International Conference on Physics in Colhsions 

(1986), ed. M. Derrick (World Scientific, 1987). 
[14] J.D. Bjorken, Phys. Rev. D 47, 101 (1992). 
[15] A.H. Mueller and W.K. Tang, Phys. Lett. B 284, 123 (1992). 
[16] O.J.P. EboU, E.M. Gregores and F. Halzen, MAD/PH-96-965 (1997). 



7 



[17] R.S. Fletcher and T. Stelzer, Phys. Rev. D 48, 5162 (1993). 
[18] E. Gotsman, E.M. Levin and U. Maor, Phys. Lett. B 3 09, 199 (1993). 
[19] E. Gotsman, E.M. Levin and U. Maor, TAUP 2485-98, |hep-ph/9804404 . 
[20] D0 Collaboration, B. Abbott et al, Phys. Lett. B 440, 189 (1998). 

[21] G. Marchesini, B.R. Webber, G. Abbiendi, I.G. Knowles, M.H. Seymour and L. Stance, Comp. Phys. Commun. 67, 465 
(1992). 

[22] G. Ingelman, P. Bruni and A. Eden, POMPYT version 2.6 

[23] A. Donnachie and P. Landshoff, Phys. Lett. B 191, 309 (1987); Nucl. Phys. B 303, 634 (1988). 
[24] CDF Collaboration, F. Abe et al, Phys. Rev. D 50, 5518 (1994). 

[25] HI Collaboration, T. Ahmed et al, Phys. Lett. B 348, 681 (1995); ZEUS Collaboration M. Derrick et al, Z. Phys. C 68, 
569 (1995). 



8 



